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In this work we derive a formalism to incorporate asymmetry and temperature effects in the 
Brown-Rho (BR) scaled lagrangian model in a mean field theory. The lagrangian density discussed 
in this work requires less parameters than the usual models with density dependent couplings. We 
also present the formalism with the inclusion of the eight lightest baryons, two lightest leptons, j3 
equilibrium and charge neutrality in order to apply the BR scaled model to the study of neutron 
- - - , stars. The results are again compared with the ones obtained from another density dependent 

' parametrization model. The role played by the rearrangement term at T=0 for nuclear or neutron 

■ star matter and at finite temperature is investigated. The BR scaled model is shown to be a good 
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tool in studies involving density dependent effective masses and in astrophysics applications. 
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^ The study of nuclear matter properties at high densities and at finite temperature has become an 
' important problem since a large variety of data, where matter is being tested at extreme conditions of 
■ density, pressure and non-zero temperature, are becoming available in the modern experimental facilities 
' I ■ which are already operational. Moreover, the constitution of the interior of neutron and protoneutron 
stars is also a problem which is receiving much attention in the recent literature. The crust of the stars. 



^ . . . 

: where density is relatively low, can be adequately described by hadronic models. The correct calculation 

^ , of the star properties as the radius and the mass depend on the accuracy of the equation of state 
(EOS) used to describe its matter. We have checked that different models, either with constant or with 



5-H ' density dependent meson-nucleon couplings present different features at subnuclear densities of nuclear 



asymmetric matter by comparing the regions of uniform unstable matter The parametrizations of 
these models generally take into account saturation properties of nuclear matter and properties of stable 
nuclei. Extensions of the models for very asymmetric nuclear matter or to finite temperatures show 
different behaviors. 

Another topic of great interest is the in-medium modification of vector meson properties. It is well 
known that the hadron masses are much larger than the sum of its constituents. One possible explanation 
for the large masses is that they may be generated dynamically Q. Spontaneous breaking of chiral 
symmetry is also related with the hadron masses. At high temperature and/or dense matter this symmetry 
is expected to be at least partially restored, which modifies the hadron masses and the decay widths 
Experimental signature of the in- medium modifications of the vector mesons have been found very 
recently |H |^ . In the in medium modifications of the uj meson were investigated in photoproduction 



experiments and its mass was found to be lowered. In [g the vector masses were verified to decrease 
by ten percent in medium. Other experimental results were also reported in the same direction Ml and 
still some experiments have been proposed to detect in-medium modifications in a near future [g. In 
order to take into account the in-medium modifications of the hadrons, density dependent relativistic 
models are certainly more useful than models with fixed parameters. In 199f Brown and Rho (BR) Q 
proposed an in-medium scaling law for the masses and coupling constants for effective chiral lagrangians. 
Our proposal here is to study hadronic properties under extreme conditions using lagrangians which 
incorporate BR scaling since this has been successfully applied to describe meson properties. As argued 
in the literature the BR scaling describes the behavior of the light mesons in extreme conditions very 
accurately. For example the enhancement of dilepton production observed in heavy ion collisions (S-f-Au) 
in the CERES collaboration and (S+W) in the HELIOS-3 is most economically and beautifully described 
by a chiral lagrangian with BR-scaled meson masses [fljlinjl. The strategy is to assume that the in-medium 
effective lagrangian has the same structure as in free space accordingly to the QCD constraints but with 
parameters which are modified in the medium. So, using this approach we obtain an effective theory with 
density dependent parameters including many-body correlations in a practical framework. It has also 
been shown that it is possible to obtain a relation between the effective parameters of chiral lagrangians 
in medium and Landau Fermi Liquid parameters . So, a link between the effective theory of QCD 
at mean field level and the many-body theory of nuclear matter is achieved. Using this reasoning the 
authors in [l^ proposed an effective Lagrangian whose parameters scale in nuclear medium according to 
the Brown-Rho (BR) scaling. 

We have already shown that one of these density dependent models, to which we refer as TW model 
[l^list . originally derived at T = can be extrapolated to finite temperatures once the thermodynamical 
consistency remains unaltered |l4l The important range of temperature which is discussed lies 

between 10 and 150 MeV since the liquid-gas phase transition takes place around 10 MeV, the phase 
transition from hadronic to quarkionic matter around 150 MeV and the relevant temperatures in the 
cooling of a protoneutron star after a supernova explosion takes place go up to approximatel y 4 MeV 
[l^ . In this work we discuss another possible application of the formalism we have derived in \lA in 
order to incorporate temperature effects in the study of lagrangians with BR scaling. 

It is well known that the same relativistic models used in order to explain data coming from heavy ion 
collisions at finite temperature, with appropriate parameter sets, also provide EOS which can be used in 
the solution of the Tolman-Oppenheimer-Volkoff differential equations ^3 for the calculation of stellar 
properties such as mass, radius and central energy density. Recent measurements of gravitational redshift 
of spectral lines provided direct constraints on the mass-to-radius ratio [isi H^. In this second case, 
however, the interpretation of the absorption features as atomic transition lines is controversial [2o[ blf. 
In recent works [23 . we have checked that there are relativistic models which can be accommodated 
within these constraints. Hence, astrophysical observations can help in the choice of appropriate models 
to describe hadronic matter. The most common relativistic model used in the description of hadronic 
matter is the non linear Walecka model [2^- When applied in nuclear astrophysics, this model is normally 
extended with the inclusion of hyperons, which are expected to appear at high densities. It was shown 
in [2^ that a low effective mass at saturation density makes the model inappropriate once hyperons are 
included. The inclusion of hyperons makes the scalar meson interaction stronger and hence the proton 
and neutron effective masses decrease more rapidly with density, acquiring a negative value. As a test of 
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the BR scaled model, we also extend it to incorporate the eight lightest baryons and enforce (3 equilibrium 
plus charge neutrality by accommodating the two lightest leptons as well. The same extension is done 
within the TW model so that both density dependent models can be compared. 

Hence, the present work is organised as follows: in section II the formalism for the BR scaled model 
and its extension to finite temperature are presented, the results are compared with the ones obtained 
with other relativist models and a discussion is included. In section III the BR scaled model is modified 
so that hyperons can be incorporated and (3 equilibrium and charge neutrality are enforced so that an 
EOS can be obtained and apphed to compact stars. In the same section the TW model is also considered 
so that the results from two density dependent models are compared. In section IV a discussion on the 
role played by the rearrangement term in different models is presented. Finally, in the last section, the 
conclusions are drawn. 



II. FORMALISM - EXTENSION TO FINITE TEMPERATURE 



In its simplest version the Lagrangian density reads 



£ = 7^ [7^ {id^ - glip) oj^ - g;{p)f ■ p^) - M*(p) + g:{p)<l>] i; 

-\p,.-r + \mfp,-fP,. (1) 

where in the notation of [3| V' is the nucleon field, the isoscalar vector field, (f) an isoscalar scalar field, 
^iiu ~ d^uju — d,yLu^, pf^ is the vector isovector field, Pf^i, = dfj,pu — dvPn — 5p(/V x Pv), t is the isospin 
operator and the masses with asterisk are BR-scaled as introduced in 

^==^ = ^ = ^ = $(P). (2) 
M mg niy nip 

The scalings of the vector coupling constants are given by 

9*s 1 9* 1 9*p 1 



(3) 



gs 1 + xp/ po' gv 1 + zp/po ' gp 1 + wp/ po ' 

where po is the nuclear saturation density. In the original papers 0, Q , g* was simply taken constant 
and equal to gs and hence, not density dependent and, with this simple choice a good description of the 
ground state was obtained. Moreover, 

Hp) = T-^' (4) 
1 + yp/po 

with y ~ 0.28 in such a way that ^{po) — 0.78 0| and z and w were taken equal or tly greater that y. 
In this work we propose another possible parametrization, based on the works |6|,|30|, where 

^p) = l-y-^, (5) 

Pa 

with y = 0.1 in such a way that a decrease of the meson masses in medium at the saturation point is 
ten percent, as found experimentally Q. It is also worth mentioning that, as far as we know this is the 
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first work where asymmetry is taken into account by the appropriate inclusion of the p meson in a BR 
scaled model. One would bear in mind that we are using the letter p for both the meson and for the total 
baryonic density. As those are common definitions in relativistic models we do not believe it may cause 
any confusion. One can see that the Lagrangian in is of the form of a Walecka-type Lagrangian and 
all the finite temperature formalism that we have developed for the density dependent hadron field theory 
(DDHFT) [Mlia can be immediately applied to these lagrangians. The thermodynamics of effective 
lagrangians with BR scaling has been studied in [i^ for zero temperature. The study of the validity of 
the BR scaling hypothesis for the non-zero temperature case is demonstrated in what follows. From 
the Euler-Lagrange equations we obtain the field equations of motion in the mean field approximation for 
infinite matter, where the meson fields are replaced by their expectation values. In this approximation, 
the expectation value of the cr, uj and p meson fields are called 00, Vq and 60 respectively. The coupled 
equations read 

rnf4>Q - g*Ps = 0, (6) 
< Vo - .g> = 0, (7) 

K^bo - ^P3 = 0, (8) 
[i^d^ - 70 {g^Vo + g;Tsbo + S^) - (M* - gj^o)] ^ = 0, (9) 

where the rearrangement term is given by 

and the scalar and baryonic densities are defined as 

Ps = (^V-), (11) 

p = (^^7^), (12) 
P3 = 2(^70r3V). (13) 

Notice that if the original parametrization for the BR-scaled model is used, g* = (/s is a constant and 
the last term of the rearrangement vanishes. In the following discussion we consider nuclear matter in 
the the mean-field approximation. Due to translational and rotational invariance the lagrangian density 
reduces to 



^-MFT = V' [iltid^' - loglVo - log*pT3bo - {M* - g*4)o)] ip 

1 *2 ,2 I 1 *2-,r2 I 1 *2,2 

"2"^^ 2"*" ^0 +2™ 



-mfq>l + -mfV^ + -mfbl (14) 



where = ±1/2 for protons and neutrons respectively. The conserved energy- momentum tensor can be 
derived in the usual fashion 24]: 



ri:^}'pT^^^i^d^i^ + g^^ 



*2 i2 -'^™*2t^2 ^^*2i2 



(15) 



Note that the rearrangement term included above and defined in eg. 1)10(1 assures the energy- momentum 
conservation, i.e., 9^T^'' = 0. From the energy-momentum tensor one easily obtains the hamiltonian 
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operator: 



where 



m^ = Ar-5:0o (17) 



P« = 2 / ^(/,+ -/,_), z-p,n, (20) 



is identified as the effective nucleon Landau mass and V is the volume of the system. Notice that the 
energy density does not carry the rearrangement term because it cancels out in a mean field approxima- 
tion: 

^ = 2$: / 70\/p^+-f (/^+ + /-) + ^^^ + ^^o^ + ^^^, (18) 

i—p,7i 

where fi^ and /i_ are the distribution functions for particles and anti-particles respectively and are 
calculated next. Following the notation in |28| . the thermodynamic potential can be written as 

Vt = £ -TS ~ ^ipPp- ^InPn, (19) 

where S is the entropy of a classical Fermi gas, T is the temperature, fXp (/in) is the proton (neutron) 
chemical potential and pp and pn are respectively the proton and neutron densities, calculated in such a 
way that p — Pp + Pn- We have 

(2^ 

where the distribution functions /i+ and /i_ for particles and anti-particles have to be derived in order 
to make the thermodynamic potential stationary for a system in equilibrium. After straightforward 
substitutions, ea. (|19|l becomes 

i—p,n 

E_/ ig. (/- >» (l^) + - M . .« {jij-) + Ml - /.-)) 
-2 E/7Srf.(/.+ -/.-)■ (21) 

i—p.n 

For a complete demonstration of the above shown expressions obtained in a Thomas- Fermi approximation 
for the non-linear Walecka model, please refer to |23|. At this point, ea. H21|) is minimized in terms of the 
distribution functions for fixed meson fields, i.e., 

if- =0 ^^J. (22) 

For the particle distribution function, the above calculation yields 



E*{p) + I]« ~p,+ g:Vo + l^&o = -T In ) ' ^^3) 
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where E*{p) = yp^ + rn]^ . A similar equation, with a sign difference is obtained for the anti-particle 
distribution function. It is important to point out that the fields 0o, Vb and 60 depend on the distribution 
function which appear in the definition of p^, p and and hence, the whole calculation is performed 
self-consistently. The effective chemical potentials are then defined as 

= M« - 5:1^0 + y6o-S^ (24) 
and the following equations for the distribution functions can be written: 

/i± = T"^ u F*! ^ ^ — ' i^P^n. (25) 

In the above calculation we have used 

and P3 ~ Pp — Pn- Within the Thomas-Fermi approach the pressure becomes 

^ = A E / , / , , (/.+ + /.-) (26) 

-^c^l + ^V,' + ^bl + i:^p. (27) 

It is worth mentioning that the thermodynamical consistency which requires the equality of the pressure 
calculated from the thermodynamical definition and from the energy-momentum tensor, discussed in 
is also obeyed by the temperature dependent Brown-Rho scaled model. 



A. Discussions on the BR-scaled model for nuclear matter 



At this point, the parameters used in the BR method have to be fixed. Through out this paper the 
nucleon mass will be M = 939 MeV, the uj meson mass my — 783 MeV and the p meson mass m„ = 763 
MeV. Three different sets are used in (27j . In what follows we use the parameter set called S3 in [23 and 
define another one which we call MA, whose bulk properties are more similar to the NL3 parameter 
set but with a larger effective mass at nuclear saturation density. While S3 is a parametrization for the 
original BR scaled model given by equation Q), MA is a parametrization for the new scaling, given in 
equation ©. In Table HI we show the S3 and MA constants and in Table Ull we display the nuclear matter 
bulk properties described by the different models used in this work. It is important to point out that the 
value for the saturation density (vide * in Table II) was not given in [iflj • For the saturation density value 
shown in Table HI the compressibility is slightly different from what is stated in (260 MeV) . This is 
probably a consequence of the fact that the authors in have not included the rearrangement term in 
their calculations. Notice that we distinguish the M* from the values. In the BR scaled models it is 
the Landau mass that should be identified with the nucleon effective mass determined by the QCD sum 
rule hxt and its value should lie in between 0.55 M and 0.75 M. 
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TABLE I: Parameter sets for the Lagrangian Q 



Set 


nis 


gs 


gv 


gp 


a; 


z 


10 


S3 


700 


5.30 


15.2 


7.97 




0.31 


0.31 


MA 


500 


7.05 


12.006 


8.761 


0.37 


0.15 


0.15 



TABLE IL Nuclear matter properties 





NL3 [M] 


TMl [33] 


GL [34] 


TW [IJ 


S3 [27] 


MA 


BjA (MeV) 


16.3 


16.3 


15.95 


16.3 


16.1 


16.3 


po (fm-^) 


0.148 


0.145 


0.145 


0.148 


0.155(*) 


0.148 


K (MeV) 


272 


281 


285 


240 


269 


258 


fsym. (MeV) 


37.4 


36.9 


36.8 


32.0 


32.0 


32.0 


M*/M 


0.60 


0.63 


0.77 


0.56 


0.78 


0.9 


ml/M 










0.68 


0.748 



In figure n we plot the dependence of the meson masses with the density for the S3 and MA parameter 
sets. This is an important consequence of this model, since the reduction of the meson masses in medium 
is an expected result [s^]. As stated in the Introduction, this behavior is related with the restoration of 
the chiral symmetry and experiments with the spectrometer HADES at GSI will also be measuring this 
effect soon. If the S3 parametrization is used, the effective masses of all mesons decrease by 22% up to 
the saturation density while if the MA parametrization is used, the decrease is forced to be just 10%, as 
found in 1^. In figure [21 we display the behavior of the coupling constants. Although both g*^s are quite 
different at subsaturation densities, they tend to achieve reasonably close values at larger densities. 5*, on 
the other hand, presents quite a similar behavior in both models and g* only changes with density within 
the MA framework. In figureOlwe show the ratios {g* /m*)'^, with i — s,v, p which are quantities always 
present in nuclear matter calculations. One can see that the ratios are very small, the ratio involving 
the uj meson being the largest in both parametrizations. One should bear in mind that, as stated earlier, 
within the original version of the BR scaled model, the scalar coupling constant does not vary with the 
density. For the MA set, on the other hand, the ratio {g*s/m*Y tends to zero at ~ 2po- 

In figure 01 we show the binding energy in terms of the baryon density for different models for T = 
and T — AQ MeV. For the sake of comparison with other models, we have chosen one model with constant 
couplings (NL3) and another one with density dependent couplings (TW). At T = 0, the TW model is 
the softest one and the NL3 the hardest, the two curves obtained with the BR scaled parametrizations 
interpolating between the other models. One can see that the temperature does not alter the softness 
(hardness) of the EOS considered. The hardest and the softest ones at T = remain so at a higher 
temperature. 

In figurej^lthe pressure versus the baryon density is displayed for symmetric nuclear matter (j/p = 0.5) 
and for very asymmetric matter (j/p — 0.1), where yp — Pp/p is the proton fraction. The isospin is a 
quantity which influences the softness (hardness) of the EOS, but one can see that the asymmetry seems 
not to affect the displayed EOSs. 

In figure El we show how the temperature affects the binding energy of the S3 and MA models for 
symmetric matter. The behavior is the same one encountered in i.e., the minimum shifts to higher 
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densities with the increase of the temperature and moves from a negative to a positive value. This seems 
to be a natural consequence of the increase in the temperature of the system. 

Another quantity of interest is the nuclear bulk symmetry energy discussed in [s^l- For symmetric 
nuclear matter at T=0 it is defined as 

= 6E*{Pp) + 8^^' ^^^^ 

with Pp = (I.Stt^p)^/'^. The value and behavior of the symmetry energy at densities larger than nuclear 
saturation density are still not well established. This quantity is important in studies involving neutron 
skins, radioactive nuclei and neutron stars. In general, relativistic and non-relativistic models give dif- 
ferent predictions for the symmetry energy. The results of this quantity for different models are also 
discussed in the present work and the values at saturation density are shown in Table ^] In figure [3 we 
plot the symmetry energy for the BR scaled models, NL3 and TW. In S3, MA and TW the symmetry 
energy at saturation density is the same. The three curves do not cross in the same point because the 
saturation densities are not the same (see Table ^J. Notice that although S3 and MA present a lower 
symmetry energy at subsaturation densities, they interpolate between NL3 and TW at larger densities. 
The value of 32 MeV that we have chosen for the symmetry energy in order to fix the gp coupling is lower 
than the ones found in most relativistic models (between 35 and 42 MeV) and approaches the values 
obtained in non-relativistic models (between 28 and 38 MeV). Notice that the choice of parameters is 
not arbitrary. They are chosen in order to reproduce the nuclear bulk properties of Table Ull Moreover, 
they also have to give the correct value of the spin-orbit splitting strength. Work in the direction of 
calculating this quantity in finite nuclei is in progress, which may require small changes in the calculated 
parameter sets. 



III. APPLICATION TO COMPACT STARS 



From the results shown in the previous section we could see that the BR scaled models show quite 
a different behavior from the NL3 and TW models at high densities. In what follows we intend to 
investigate which are the consequences of using the BR scaled model in the description of neutron star 
matter. The behavior of the EOS at high densities is responsible for the determination of the maximum 
mass of the star. In order to apply the BR scaled density dependent model to compact stellar objects, 
it is important to allow for the inclusion of the eight lightest baryons (nucleons. A, S°, E^,S~ and S°) 
as well as the two lightest leptons (e~ and /i). The baryons have to be considered since their masses are 
such that their presence is already possible at the neutron stars high densities. The leptons, on the other 
hand, play a decisive role in ensuring charge neutrality and f3 equilibrium. As seen in equation (|24|l the 
rearrangement term which appears due to the density dependent couplings alters the chemical potentials 
of the particles in the system and hence, the /3 equilibrium conditions are somewhat different as compared 
with the usual relativistic models. In what follows we show the formalism developed for the BR scaled 
model in neutron stars and also for the TW model, so that two density dependent parametrization models 
can be compared. 
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A. Considering (i equilibrium within the BR-scaled model 

For our purposes of testing the BR scaled model at high densities, which are present in neutron stars, we 
shall restrict ourselves to the T — case. Of course, the extension to finite temperature is trivial and can 
be done following the steps of section II. Actually in calculations involving protoneutron stars or stars 
with fixed entropy and trapped neutrinos, the extension has to be done. Equation is then modified 
in order to accommodate these new particles 

B 

-\p^. ■ + \^fp^ • + E (^^'-^^ - "^') ■ (29) 

I 

where the meson field operators represent the same mesons as in eq.lQ), '0b now represents each of 
the eight baryons, I describes the two leptons whose masses are respectively — 0.511 MeV and 
= 106.55 MeV and the masses with asterisk are again BR-scaled as in eq. 0, the 8 baryons of the 
octet also obeying the same scaling law, i.e.. 

The baryon meson couphngs are defined as g*^ = x^b Qs, 9vB = ^-"B 9^, 9*pB = ^pB 9*p and XsB-, x^b 
and XpB are equal to 1 for the nucleons and may have different values for the hyperons. Again the meson 
fields are obtained in the same way as in section II and they now read: 

mf(l^o-^9*sBPsB = 0, (31) 

B 

mfVo-Y^glBPB = 0, (32) 

B 

^f^0-^9*pB'r3BPB 0, (33) 
B 

[^rO^ - 70 (gM + 9;b r3B bo + ^^br) - iM*B - 9:bM] ^ - 0, (34) 
where the term S^^^ is now changed and is given by 

B B 
B 

T^B is the isospin projection of each baryon and the scalar and baryonic densities are 

PsB = ^/ P'dp^, (36) 
PB ^ fi, (37) 
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with 

^LB = M*g - g*B4>o, (38) 



E*B = ^P^+mls^ (39) 

VB = MB - .g^ijV^o - g*pB T3B bo-^oBR (4^^) 

kU+^lB^- (41) 

The equation of state, necessary for the description of the stellar matter can now be obtained. The 
energy density and the pressure density are given respectively by 

^-^E/ P'dpE^B + '^^^+'^V.' + ^bl + ^Y. P'dpEi (42) 
^ B -IQ ^ ^ ^ ^ I 

and 

► 2 m*^ T r^Fi ^4, 



(43) 



where Ei = ^/p"^ + rnf , pi = and fci?; = fif — rnf, /i/ being the chemical potential of the lepton I. 
Notice that, as far as the leptons do not exchange mesons with the baryons nor with themselves, they 
were introduced as free Fermi gases. The weak interaction between leptons and hadrons is taken into 
account through the constraints of charge neutrality and /3 equilibrium given respectively by: 

E^sPB + ^gfPi^O, (44) 

B I 

where Qb is the electric charge of baryon B, qf is the electric charge of lepton I and 

^J■B = ^J■n - q^BlJ-e- (45) 



B. Considering (5 equilibrium within the TW model 



In order to make a comparison with the density dependent BR scaled model, we next make some 
considerations about the TW model (l^ , originally derived at T = and which has also been extrapolated 
to finite temperatures 0, . In what follows we write the most important formulae for the TW model 
once the lightest baryon octet and the lightest leptons are included and charge neutrality and [3 equilibrium 
are enforced. In reference [s^l a similar approach was developed and two different models were discussed. 
In the first of them the couplings depend on the total baryonic density, as done in the BR scaled approach 
shown in the last subsection and also next in the present work. In the second model the couplings depend 
only on the proton plus neutron densities. The authors of reference js^ showed that an examination of 
the neutron star properties favored the first model. Notice that we next redefine many of the previously 
defined quantities, as the baryon effective mass, baryon chemical potentials, etc. The new equations 
should not be mixed up with the equations given in the previous subsections although we have kept the 
same notation. 
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We start from the lagrangian density 

^ = ^^s[7M(«5''-r„BCc;^-rpsf.p^)-(Af-r,B</))]VB 

B 

■ + \mlp^ ■f^ + Y.^i (*7m5'^ - mi) (46) 
I 

with all the definitions for the fields given after eq. p9|l still holding. Yib and are respectively the 
couplings of the mesons i — s,v, p with the hyperons and their masses. In this model the set of constants 
is defined by Tsb — x^b Ts, T^b — x^b Tu, ^ pB = Xps Tp and as in subsection III A, Xsb, x^b and 
XpB are equal to 1 for the nucleons and can acquire different values for the hyperons. Fg, and Tp 
are the nucleon-meson coupling constants which are adjusted in order to reproduce some of the nuclear 
matter bulk properties, using the following parametrization: 

^i{p) ^^^{Po)f^{x), i = S,V (47) 

with 

^^^")="^l + c.(. + dO- ^''^ 

where x = p/po and 

Tp{p) - Tpipo) exp[-ap(x - 1)], (49) 

with the values of the parameters rrij, Tj, aj, bi, Ci and di, j — s,v, p given in |l2| . The nuclcon, uj and 
p meson masses are taken as in the BR scaled model. The scalar meson mass iris is 500 McV. Other 
possibilities for these parameters are also found in the literature j^^. The meson and baryon coupled 
equations for the fields read: 

ml(f>o - ^ TsB PsB = 0, (50) 

B 

mlVo-^VyBPB = 0, (51) 

B 

TOp&o - X! '^^^ ^ ^^'^^ 

B 

[il^'dp, ~ 70 {T.bVo + TpB T3B bo + ^^tw) " ^b] = 0, (53) 

where the term S^^,^ is given by 

\-R \dTyB , dVpB , dVsB , 
^0 TW = ~Q^P^^° + ~dp~'^'^^ dfr^""^ 

and the scalar and baryonic densities are defined as 



(54) 



1 r M* 

P^B = ^ j P'dp^, (55) 
PB = (56) 
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with 

M*B = Mb^TsbcIjo, (57) 



E*s = ^p2 + M*2. (58) 

Notice that the rearrangement term shown in equation H54I) is the same one shown in equation (18) of 
[s^ . once the delta and the strange mesons are excluded from their calculation. The effective chemical 
potentials are then defined as 



i/B = MB - TvbVo - TpB T3B bo - S^Vw ^ \J ^Ib + ^^B ■ (59) 

The conditions for P equilibrium and charge neutrality are again the same ones as in subsection III A, 
given by eas. H44|l and H45|) . The final expressions for the energy density and pressure become respectively: 



~ g Jo ^ ^ ^ ~ J Jo 

and 

rkpB ™2 ^2 ^2 -, „kFi 



(60) 



1 ^ [^-- p^dp mi . mt . „ 1 ^ f^" p^dp 



C. Discussions on the compact star properties 

Although NL3 0| and TMl 0| are the most common parametrizations of the NLWM for nuclear 
matter and finite nuclei studies, it is well known that they are just adequate for the description of neutron 
star properties if only protons, neutrons and leptons are considered as possible constituents 25]. The 
inclusion of hyperons softens the EOS, but also makes the baryon effective masses decrease too fast and 
a good convergence can only be obtained at relatively low densities. For this reason, whenever hyperons 
are considered in the present work, we shall make comparisons with the GL js^ parametrization of the 
NLWM, where the above mentioned problem does not exist. It is our aim also to verify whether this 
problem is present in the TW and BR scaled model. 

At this point the meson-hyperon couplings have to be fixed. Several possibilities are discussed in the 
literature According to QjE^ the hyperon couplings constrained by the binding of the A hyperon 

in nuclear matter, hypernuclear levels and neutron star masses yields XsB = 0.7 and x^b = XpB ~ 0.783 
and the couplings to the S and S are equal to those of the A hyperon. Another possibility is to take 
XsB = XvB — XpB — ^2/3 as in [SEEE^I- This choice is based ^on quark cou nting arguments. The 
universal coupling, with Xsb = x^b = XpB = 1 has also been used ■40j. From |23i l4l| it can be verified 
that the compact star properties depend on the choice of these parameters. As our aim in the present work 
is to compare results obtained from different models and the correct choice is still not well established, 
we have used the simple universal coupling in what follows. 

In figure IHl we show the particle population for the NLWM with the GL parametrization, for the TW 
and the BR scaled model either with nucleons only or with the 8 baryons. The particle fraction is defined 
as Yi — Pi/ i — ^ baryons and 2 leptons. We have again chosen two possibilities for the BR scaled 
parameters, namely S3 and MA. If just protons and neutrons are included, the TW model presents a 
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slight decrease in protons and consequently slight increase in neutrons and the MA model shows the 
opposite behavior as compared with the GL parametrization. If hyperons are considered the results are 
all quite different, as a consequence of the different EOS shown in figure IHl One can clearly see that 
all surviving particles tend to the same amount in the BR scaled model, probably a consequence of the 
enforced scaling law, a feature which happens earlier within S3 than with MA. From figure |^ one can 
see that in both cases, the TW EOS is the softest one and the MA the hardest. S3 and GL interpolate 
between the other two EOS. For the present choice of parameters, the TW model also stops converging 
at a too low density for astrophysical studies. The consequences of this fact will be discussed next. 

Once the EOS are obtained, we solve the Tolman-Oppenheimer-Volkoff equations in order to obtain 
the stellar properties. They read 

dP ^ G [£ + P] [M + Anr^P] 

dr r {r-2GM) ' ^ ' 

^ = inr'S, (63) 
ar 

with G as the gravitational constant and M(r) as the enclosed gravitational mass. We have used c — 1. 
Given an EOS, these equations can be integrated from the origin as an initial value problem for a given 
choice of the central energy density, (^o)- The value of r {= R), where the pressure vanishes defines the 
surface of the star. In Table IIIII we display the results for the stars with the maximum gravitational 
mass, the maximum baryonic mass, their radii and central energy density for each of the EOS discussed 
in the present work. One can see that, if only nuclcons arc considered, the TW model presents the lowest 
maximum mass and the smallest radius with a consequent very large central energy density. If hyperons 
are included the obtained result for the TW model is just shown for completeness because it is unrealistic 
once the maximum mass was not achieved since the program failed to converge at high densities. The 
results obtained for the maximum masses with the BR scaled models are somewhat larger than with the 
GL model, but still in the expected range of values. Different results can be obtained with a different 
choice of the x^b, x^b and XpB constants. Once hyperons are included, the EOSs always become softer 
with a consequent lower value for the maximum stellar masses and radii and larger central energy density. 



TABLE III; Hadronic star properties for the EOSs described in the text 



type hadron model Mm.ax{MQ) Mb maxiM^) i? (Km) eo (fm"") 

np GL 2.40 2.89 12.19 5.43 

np+hyperons GL 2.18 2.56 11.35 6.34 

np TW 2.08 2.46 10.62 7.20 

np+hyperons(*) TW 1.89 2.24 9.46 8.44 

np S3 2.88 3.57 12.81 4.57 

np+hyperons S3 2.65 3.36 11.33 5.60 

np MA 2.86 3.59 11.79 5.29 

np+hyperons MA 2.76 3.49 11.00 5.93 
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IV. THE IMPORTANCE OF THE REARRANGEMENT TERMS 

In what follows we concentrate on the role played by the rearrangement term in the different models. 
This term is very important and influences all properties of nuclear and neutron star matter since it 
changes the effective chemical potentials given in equations (I24f) . (|41l) and (|59|l . 

In figure IIUI we show the rearrangement term which arises in various situations. In figure llUb . it is 
shown the rearrangement terms which appear in symmetric nuclear matter (without the imposition of 
f3 stability) for the TW and the two parametrizations of the BR scaled model within the density range 
considered in nuclear matter studies. The S3 parametrization produces a term which is much more 
negative than the TW model. The MA parametrization shows a rearrangement term which decreases 
even further. So, in this case the rearrangement term is more attractive within the BR-scaled models. 
Notice that in the S3 parametrization, if y and z were chosen as having the same value, only the scalar 
meson would contribute to the rearrangement term. 

Please, notice that the scales of the figures mentioned next are all different. In figures [TUb and 1 10b we 
show the influence of the temperature on the rearrangement term of the two BR scaled models. As the 
temperature affects very little the term in the MA parametrization, the same does not happen if the S3 
parametrization is chosen. Moreover, the rearrangement term increases slightly as temperature increases 
within the MA and decreases quite a lot with the increase of temperature with the S3 choice. 

In figure llOL l. the rearrangement term of the TW model is shown for symmetric nuclear matter, for 
a very asymmetric nuclear matter, with a proton fraction yp = 0.1 and for the equation of state where 
charge neutrality and f3 stability is required either with only protons and neutrons or with hyperons as 
well. The same is shown in figureslTUfe andllOt' for S3 and MA. While in the TW model the rearrangement 
term decreases with the asymmetry of the system and increases when the conditions of f3 equilibrium and 
charge neutrality are enforced, the influence of the hyperons being very small, in the S3 all curves are 
very similar and in the MA the asymmetry almost does not interfere in the rearrangement term and the 
hyperons again play no role, but charge neutrality and /3 equilibrium conditions modify the curve quite 
drastically. 

In figures lint' and II Oh one can see a comparison between the rearrangement term arising from the TW 
model and the ones obtained with the BR scaled models respectively for the case when only nuclcons are 
considered and when hyperons are also included in the system. We have now opted to show a much wider 
density range so that the differences can be clearly seen. While the inclusion of hyperons makes the TW 
model rearrangement term increase slightly at densities of the order of ~ 2 fm~'^ as compared with the 
case when only nucleons are considered, at about the same density the MA rearrangement term starts to 
decrease. The S3 term also decreases more rapidly if hyperons are considered, but it starts at a density 
of the order of ~ 1 fm^'^. At these very large densities other important features as the deconfincment to 
the quark matter is already present. 

V. CONCLUSIONS 

In the present work we have derived a formalism to incorporate temperature effects in the BR scaled 
model to make it useful in future heavy ion collision studies. We have also investigated the possibility of 
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applying it to nuclear astrophysics by enfor cing; charge neutrality and /3 equilibrium. We have compared 
our results with the more standard NL3 |29j version of the NLWM in nuclear matter and with the 
GL [s^ parametrization of the NLWM in neutron star matter. In both cases, the BR scaled model 
was also compared with another density dependent model, the TW Q. It is worth pointing out that 
density dependent models are alternative approaches to describe hadronic matter without the usual 
non-linear terms necessary in the NLWM. Although in the low baryonic density regime all models and 
parametrizations used in this paper have EOSs with similar behavior, different scenarios show up when 
the density increases, specially, when the baryon octet is taken into account. The lagrangians with BR 
scaling are much simpler than the other model with density dependent couplings (TW) and they provide 
effective meson masses which decrease with the increase of the density, behavior which has already been 
confirmed by experiments Q and which is not present in all models where the meson masses are 
held fixed. We have also shown that the BR scaled model is a very good tool in describing neutron 
star properties. Of course a more systematic study can be done by including the delta E^l and the 
strange ^3 mesons in the lagrangian density. Other possible choices for the meson- hyperon couplings 
should also be considered. More realistic hybrid stars with a deconfinement to the quark phase can now 
be obtained within the BR scaled model for the hadron phase. When the new EOSs are built and the 
stellar properties are obtained, the mass to radius constraints Q Q can be used as a probe to the 
new formalism. It is also important to say that the rearrangement term plays a central role in density 
dependent models and many of the system properties depend on its strength.lt is crucial to assure the 
energy-momentum conservation and the thermodynamical consistency for density dependent models. In 
contrast to the pressure and compressibility which depends explicitly on the rearrangement term, in the 
energy density it cancels out. However, it still contributes to the binding energy through the chemical 
potentials. In figure [Till it was shown that the BR parametrizations yields more attractive rearrangement 
term as compared with the TW model. 

We have checked that other possibilities for the parametrization of the BR scaled model given by 

£ ^ 1 + xip/po 9l ^ 1 + x-ipjpQ ^ l + x5p/po ^g^^ 
gs l + x2p/po' gv l + X4p/po' gp l + a^ep/po' 

also work. Nevertheless, one of the advantages of using the BR scaled models instead of the TW model 
is that they contain a much smaller number of parameters. With the choice given in (|54|l much of the 
beauty of the model would be lost, but it remains as a possible alternative. Finally, it is important to 
stress that a test to finite nuclei in order to obtain the correct value for the spin-orbit splitting strength 
is still necessary with the BR parameter sets used in this work. 
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FIG. 1: Meson masses m* ,i — s,v,p versus the baryonic density. 
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FIG. 4: Binding energy in terms of the baryon density for different models at T = (top figure) and T = 40 MeV 
(bottom figure). 




FIG. 7: Symmetry energy in terms of the baryon density for different models at T 
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FIG. 8: Particle fractions for a) the non linear Walecka model with the GL parametrization with nucleons only, 
b) with 8 baryons, c) TW model with nucleons only, d) with 8 baryons, o) BR model with parameter sets S3 and 
nucleons only f) S3 and 8 baxyons, g)MA and nucleons only, h) MA and 8 baryons. 




FIG. 9: EOS for the models discussed in the present work with nucleons only (top figure) and 8 baryons (bottom 
figure) . 




FIG. 10: Rearrangement terms calculated in different situations: a)difFerent models in symmetric nuclear matter, 
b) S3 and symmetric nuclear matter for different temperatures, c) MA and symmetric nuclear matter for different 
temperatures, d) TW, e) S3, f) MA, g) nucleons only, h) Sbaryons. 



